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Abstract. The problem of the self-interaction of a quasi-rigid classical particle with 
an arbitrary spherically symmetric charge distribution is completely solved up to the 
first order in the acceleration. No ad hoc assumptions are made. The relativistic 
equations of conservation of energy and momentum in a continuous medium are used. 
The electromagnetic fields are calculated in the reference frame of instantaneous rest 
using the Coulomb gauge; in this way the troublesome power expansion is avoided. 
Most of the puzzles that this problem has aroused are due to the inertia of the 
negative pressure that equilibrates the electrostatic repulsion inside the particle. The 
effective mass of this pressure is — U c /(3c 2 ), where U c is the electrostatic energy. When 
the pressure mass is taken into account the dressed mass m turns out to be the 
bare mass plus the electrostatic mass to = too + U c /c 2 . It is shown that a proper 
mechanical behaviour requires that too > U e /(3c 2 ). This condition poses a lower 
bound on the radius that a particle of a given bare mass and charge may have. The 
violation of this condition is the reason why the Lorentz-Abraham-Dirac formula for 
the radiation reaction of a point charge predicts unphysical motions that run away or 
violate causality. Provided the mass condition is met the solutions of the exact equation 
of motion never run away and conform to causality and conservation of energy and 
momentum. When the radius is much smaller than the wave-length of the radiated 
fields, but the mass condition is still met, the exact expression reduces to the formula 
that Rohrlich [ Phys. Lett. A 303 307 (2002)] has advocated for the radiation reaction 
of a quasi-point charge. 
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1. Introduction 

Perhaps the last unsolved problem of classical physics is the self-interaction of charged 
particles. We refer to a classical particle with a spherically symmetric charge density 
moving under the influence of an external force F ex and the force F em due to the 
electromagnetic fields that it itself generates. This last self-force can be expressed as 
the sum of the negative time derivative of the momentum of the local electromagnetic 
fields that surround the particle (Fl) plus the radiation reaction (Fr). Fairly complete 
relations of more than a century of history can be found in the monograph of Yaghjian 
PP, in the review of Rohrlich [2j and in the chapter 16 of the textbook of Jackson [Hj. 
The published solutions of this problem are filled with puzzles. Let us consider the most 
notorious ones. We will use Gauss units. 

1) Spontaneous self- acceleration. The Lorentz-Abraham-Dirac (LAD) formula for 
the radiation reaction of a point charge |H E], which in the non-relativistic limit 
reduces to -Fr = |g 2 a/c 3 , has solutions in which the particle, in the absence of an 
external force, spontaneously accelerates itself approaching speed c (run-away solutions). 

2) Violation of causality. There are also solutions of the LAD equation that do not 
run away, but where the acceleration depends on values that the applied force will take 
in the future. 

3) Power unbalance. The electrostatic field (i. e. the electric field that is present 
when the particle is not accelerated) does not contribute to the net force on the particle, 
because the integral over the volume of the corresponding force density vanishes. It, 
nevertheless, contributes to the power, because of the Lorentz contraction. As a result 
the total power of electromagnetic forces differs from the product F em - v. This was 
noticed by Abraham in 1904 In 1905 Poincare [HUH] showed that the existence of a 
pressure was required in order to balance the electrostatic repulsion and thus assure the 
stability of the particle. The work of the pressure as the volume of the particle changes, 
was supposed to fix the power unbalance. But there is a problem, in Ref. jH] and in 
the introduction of Ref. [0] (p. 20) Poincare refers to an external pressure, whose 
work is indeed proportional to the change of volume, but the pressure that is actually 
calculated (p. 63 of Ref. [5]) is an internal negative pressure, whose work over the whole 
particle is obviously zero. It is worth remarking that the results of Lorentz, Abraham 
and Poincare on this subject conform with relativity in spite of being previous to the 
1905 Einstein paper, because they assume an ad hoc Lorentz contraction. 

4) The 4/3 problem. The electrostatic energy U e should contribute to the mass of the 
particle, so if mo is the bare mass then the dressed mass must be m = mo + U c /c 2 . But 
if one computes the momentum of the electromagnetic field that surrounds a spherical 
charge distribution moving with velocity v one obtains ^U^v/c 2 . This value agrees with 
that part of the self- force that is proportional to the acceleration (-FL). To make matters 
more confusing the energy of such field is U e j(l + \{. v l c ) 2 )i which has an additional 
|(w/c) 2 term. 

Schwinger has shown, for a spherically symmetrical charge distribution moving 
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uniformly, that there are tensors which depend on the fields, that can be added to the 
electromagnetic energy-momentum-stress tensor T^ v in order to obtain a divergence-less 
total tensor. He proposed that such "corrected" tensors should be used for calculating 
the four-momentum of the fields. He found two possible tensors, corresponding to 
electromagnetic masses U e /c 2 and |£7 e /c 2 . Rohrlich [2] has remarked that any value of 
the electromagnetic mass in between is also possible. 

Instead Singal JT] has shown that the conventional theory is consistent when 
a proper account is taken of all the energy and momentum associated with the 
electromagnetic phenomenon in the system. 

5) Apparent violation of energy and momentum conservation. The Larmor ^2] 
formula for the radiated power by a point charge in the instantaneous rest frame is 
-Prad = f<7 2 a 2 /c 3 . But the LAD formula predicts that if the particle is moving with 
constant acceleration the reaction force is zero and therefore also its power. From 
where is the radiated power coming? A similar situation occurs with the momentum. 
In the rest frame the rate of radiated momentum vanishes. Where does the impulse 
of the radiation reaction force go? One should expect that the electromagnetic field 
surrounding the particle acts as a reservoir of energy and momentum. 

The Classical Electromagnetism of a continuous current density is a perfectly 
well behaved theory that conserves energy and momentum and conforms to causality. 
Solutions that run away or violate causality simply cannot appear. Various are the 
sources of the problems. 

The run-away solutions and the violations of causality should be attributed to the 
R — > limit. We will show for any spherical charge distribution, that in order to avoid 
run-away solutions the total effective mass of the particle m must be greater than the 
mass of the electromagnetic fields surrounding it, i. e. ^U e /c 2 . The reason is simple: 
m — ^U e /c 2 is the effective mass of matter, if it were negative the acceleration would be 
opposite to the force. As U e > \q 2 / R there is an absolute minimum for the radius of a 
classical particle of mass m and charge q, R > \q 2 1 (mc 2 ). Therefore the R — > limit is 
unphysical, because R = is interior to an unphysical region. 

The violation of causality is due to the expansion in powers of R that is used to 
evaluate the electromagnetic fields and the self-force. As the fields are linear functionals 
of the retarded currents j M , so must also be -F cm - As this self-force vanishes if the particle 
is not accelerated, it should be a functional of the retarded acceleration a(t — \x — y\/c), 
where x and y are two points of the particle. We will find such functional for an 
arbitrary charge density. So, for the exact solution of an extended particle there are no 
pre-accelerations or violations of causality. Of course if we make an expansion in powers 
of \x — y\ and we take any finite number of terms we get an instantaneous result which 
violates causality. If we take the R — > limit we are left with the R~ l and the R° terms, 
the first one diverges and corresponds to the local fields reaction Fi, while the second 
is the LAD result for the radiation reaction as said above such limit is unphysical. 
In the calculation of Dirac |Hj the violation of causality is explicit, because the advance 
Green's function is used. 
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The exact self- force of a sphere of radius R with a uniformly charged surface has 
been calculated by Sommerfeld [13] . Page [E|, Caldirola [T3] and Yaghjian [1J. Their 
result confirm what we have been discussing [21 EE]. Our general result yields the 
Sommerfeld-Page- Caldirola- Yaghjian formula for the charged spherical shell. 

The power unbalance and the 4/3 problems have a common origin. They are due to 
having overlooked a purely relativistic effect: pressure has inertia. In our calculation we 
use the relativistic momentum-energy conservation formula for a continuous medium. 
Here the momentum-energy tensor of matter, the stress tensor and the force density 
appear. We have to remark that this stress tensor is not some ad hoc addition, it is just 
the stress tensor that any continuous medium should have. At rest the stress tensor is 
determined by the electrostatic repulsion and is purely spatial, but when the particle 
is moving the temporal components of the tensor do not vanish. Those components 
contribute to the momentum density and to the energy density of the particle. By 
integration over the volume one obtains the equation of motion and the energy equation. 
The pressure contribution to the mass of the particle is — \U e /c 2 which is just what is 
needed to fix the 4/3 problem. The power unbalance also disappears. 

Finally the apparent violation of energy and momentum conservations is understood 
by considering all the contributions to the energy and momentum of the electromagnetic 
field that surrounds the particle. When the particle is not accelerated there is an 
electromagnetic field around the particle, the local field, which contributes to the energy 
and momentum of the particle. Call and its energy density and Poynting vector. 
In addition to this, when the particle is accelerated there is also the radiated field with 
energy density «r and Poynting vector Sr. But the energy density and the Poynting 
vector are quadratic functions of the fields, so the total energy density and Poynting 
vector have additional cross terms, S = + Sr + Sq and « = m l + « R + « c These 
cross terms act as energy and momentum reservoirs. 

The point charges (R = 0) are inconsistent with classical electrodynamics, but 
phenomenologically a point charge is a body with a radius that is negligible as compared 
to the other relevant distances, in particular to the wavelength of the fields with which it 
interacts. We will call such particles physical point charges, while we call mathematical 
point charges those with R = 0. We will show that the acceleration can be expressed as 
a time integral of the retarded applied force. When the physical point charge conditions 
are fulfilled the expression reduces to an equation similar to LAD motion equation, but 
replacing a by F ex /m, exactly as was recently proposed by Rohrlich [T7j . 

2. Quasi-rigid motion 

A perfectly rigid body is inconsistent with relativity because it implies instantaneous 
propagation of signals. In general if a body is accelerated it is stressed and deformed. 
How this happens depends on the particular characteristics of the body and the applied 
force. If the external force is more o less uniform through the body and the acceleration 
is small and does not change too rapidly the body can, as seen in the instantaneous rest 
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frame, continuously maintain its shape suffering a negligible deformation compared with 
particle radius. The bound on the acceleration a is easily determined. The speed change 
during the time that a signal takes to traverse the particle should be negligible with 
respect to the speed of sound (speed of signals in the particle). So, {2R/v s )a <C v 8 . For 
a particle of maximum stiffness; v s — c and so a < c 2 / (2R), which is an immense bound 
even for macroscopic particles. The bound of the time derivative of the acceleration is 
obtained by comparison of the relative rate of change of the acceleration with the inverse 
of the transit time, a/a <C v s /(2R). Again assuming maximum stiffness, a/a <C c/(2R). 
If both conditions are satisfied the motion of the particle is essentially rigid and the 
internal degrees of freedom of the particle are not excited. We will assume that this is 
the case. We are explicitly excluding effects produced by deformation of the particle, 
such as polarisation. 

Consider a spherically symmetric particle. Let Xo(t) be the position of the centre 
of the particle in the laboratory reference frame. Let v (t) and a (t) be the velocity and 
acceleration of that point. The instantaneous rest frame would be defined by a Lorentz 
transformation of velocity v and "parallel" axes. Let y be the position relative to the 
centre of an element of the particle in the rest frame. We will assume that the motion 
of each element of the particle is given by the following expression 

x{y,t) = x (t) + y + (% x - l)(y ■ v )v , (1) 

were 70 = (1 - (v /c) 2 )-2. This expression is equivalent to the Lorentz transformation 
of the laboratory relative position x(y,t) — xo(t) to the rest frame position y. If the 
particle is not accelerated (0) reproduces correctly the Lorentz contraction. 

If the particle is accelerated the velocity is different in different points of the particle 

v (y, t) = MM (2) 

= v (t) + ^-(h ■ «„)«„ + ( 7o -' - ■ «„)«„], (3) 

but in the rest reference frame all the points are at rest v(y,t) = 0. The accelerations in 
different points differ even in the rest frame, but the difference in that case is negligible 
a(y,t) = (1 -y a /c 2 )a . 

It is convenient to define Sv as the difference between the velocity of a point and 
the velocity of the centre 

5v = v(y,t)-v (t). (4) 



3. Extended particle equation 

The dynamics of an extended particle is determined by the momentum-energy tensor 
O^, the stress tensor P^ v and the force density f 1 with the equation [TH] 

v At (e^ + p^) = p. (5) 

For a spherically symmetric particle the proper density of rest mass fi(y), which 
is a four-scalar, is a function of the radial y = \y\ coordinate in the rest frame. The 
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integral of fi(y) over the proper volume is the bare rest mass of the particle m . The 
tensor is given by 

= /juW, (6) 

where w M is the four- velocity (V = 71/, u° = 7c). 

The stress tensor in the rest frame is symmetric and purely spatial, P 0u = P u0 = 0. 



We are using the metric tensor g^ v that has a positive trace. In general P 00 = v(VjP % i /c 2 
and P°* = P' l ° = P^Vj/c. At rest, because of spherical symmetry, the stress tensor must 
be of the form P = P r rr + P t (I — rr), where I is the identity tensor and r is the radial 
unit vector. Also at rest the force density / is the electrostatic repulsion and the (JHJ) 
reduces to V • P = /. This equation determines P, but the relation between the radial 
stress P r and the transverse stress P t depends on the particular elasticity properties of 
the particle. Nevertheless, the final result after integrating over the volume, should be 
independent of the particular elasticity model. Therefore, for the sake of simplicity we 
will assume a Pascalian (isotropic) pressure; that is, in the rest frame P u = P5ij. The 
pressure P is a four-scalar. Then in a generic reference frame 

pv» = P(g^ + u^u u /c 2 ). (7) 

The momentum and energy equations are then 

^[(/x + P/c 2 ) 7 2 v] + V ■ [(// + P/c 2 )j 2 vv] + VP = / (8) 

and 

J^K/ic 2 + P) 7 2 - P] + V ■ [(/ic 2 + P) 7 2 ^] = f v. (9) 

It is interesting to compare these equations with the non-relativistic ones. The only 
pressure terms that appear in the non-relativistic equations are the gradient in (JHJ) and 
the pressure that is inside the divergence of Q . Note that the Pj c 2 terms in (jHJ) do not 
disappear in the limit v/c — > 0. This means that, in this case, the usual non-relativistic 
limit, v/c — > 0, does not yield the actual non-relativistic result. In other words, the 
inertia of pressure is a purely relativistic phenomenon. 

Neglecting the inertia of pressure is at the origin of some of the discussed puzzles. 
The P/c 2 terms that add up to \x are related to the 4/3 problem. The additional — P 
term in the time derivative of JHJ) is related to the power unbalance. 

In order to obtain the equations of motion for the particle we integrate over the 
whole volume. As \i and P vanish outside the particle the divergence terms disappear. 
That is the reason why, in the Newtonian mechanics, the pressure has no effect on the 
motion of the whole particle. Let us consider the first term of the left hand side of (JHJ), 

d 3 x -[(/! + P/c 2 ) 7 2 ^] = -Jd 3 x(fi + P/c 2 h 2 (v + 5v) (10) 
d 



-[ 7o 2 / d 3 x(fi + P/c 2 )(v + 6v)] + O(a 2 ) (11) 
[ 7o 2 / d 3 x(fi + P/c 2 )v ] + O(a 2 ) (12) 



d 
dt 
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= ^[7o / d 3 y Oi + P/c> ] + 0(a 2 ). (13) 

In (JTUJ) we have used (jlj) In (JTTJ) we have replaced 7 by 70, 7 = 70 + 0(a). In (fTTj) 
the term of 5?; vanishes because 5v is odd in y. In (llHj) we have transformed the integral 
to the rest frame, which gives a factor 7^ .At the end we are left with 

^[(m + mp)7o-y ] = F (14) 
where we have defined the pressure mass m P and the total force F as 
mp = c~ 2 

and 

F 

Similarly, 

d 



/ d 3 yP (15) 
Jd 3 xf. (16) 



[(m + m P )c 2 7o - mpc 2 ^ 1 ] = F ■ v + J d 3 x f ■ 5v. (17) 



dt 

Because of the following identity 
2 d7o d7 v 

if mp is constant, it must be that 
2 d 7o _1 



— mpc 



(1/ Jd 3 xf-5v. (19) 

The pressure has two effects: First it contributes to the effective mass of the particle. 
Second, as long as (fT9"j) is verified, it cures the power unbalance. The pressure mass mp 
could be negative, as it is in the charged particle case. A proper mechanical behaviour 
requires a positive effective mass m + mp > 0, otherwise the acceleration would be 
opposite to the force. 

In the rest of the paper there will be no confusion between the velocity of a point 
in the particle and the velocity of its centre, so in the following we will drop the zero 
subscript of v , a and 70. 

4. Stresses in the charged particle 

Consider a particle with a spherically symmetric charge density, that in the rest frame 
is p(y) = qg{y), where q is the charge. The fraction Q(r) of charge contained inside a 
sphere of radius r is defined as 

Q(r) = 4vr [ dyy 2 g(y); (20) 

of course 







^ = Wg{y) (21) 
dy 
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and g(y) is normalised to 1, 



J d 3 yg(y) = BmQ(j/) = 1. 



(22) 



The convergence of the integral requires 
Urn y 3 g(y) = 0. 

It is also required that the electrostatic energy be finite, this implies 



(23) 



\im y 5/2 g(y) = 



(24) 



and 



lim J /- 1 / 2 Q(y) = 0. 



(25) 



The force density is the sum of an external force density and the contribution 
of its own electromagnetic fields, / = / ex + f cm , with f em = pE + c -1 j x B. The 
electromagnetic field F^ v has two components, a local field F£" which decays as r~ 2 
and the radiation field F£ u which decays as r _1 and is proportional to the acceleration. 
The self-force / cm has therefore four terms corresponding to the electric and magnetic 
forces of the local and radiation fields. The local electric field depends on acceleration 
Ei,(a). We will define the electrostatic force density as / es = pE^{0) ; this is the 
largest term, by definition it is independent of acceleration, but by symmetry it does 
not contribute to the net force. It does contribute to the power. The remaining local 
electric force density p(E L (a) — E L (0)) and the other three terms will be included in 
the field reaction force density /f r , so f em = f es + /f r . The force of the local magnetic 
field, of course does not contribute to the power, but it does give a contribution to the 
total force which is of first order in the acceleration. 

The momentum and energy of the local field are bound to the particle and 
contribute to its total energy and momentum. When the particle moves with constant 
velocity the local field is easily obtained with a Lorentz transformation of the electric 
field at rest. When the particle is accelerated the local field is modified far away from 
the particle. So the local field momentum and energy are only slightly dependent on 
the acceleration. The zeroth order can be calculated using the constant velocity fields. 

At rest the electric field is 



The electrostatic energy U c can be evaluated with any of the following equivalent 
expressions 




(28) 



(27) 



(29) 
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At rest the (JSJ) reduces to VP = / es , from where the pressure can be obtained. In 
general it also depends on / ex and on / fr , but as f\ r <C / es and as it is supposed that 
the external force does not deform the particle, those dependencies will be neglected. 
In the frame of instantaneous rest the pressure, up to zeroth order in the acceleration, 
is then 



P(r) = -<rl dy 9iy)Q 2 iy) +0( fl ). (30) 



m P = ~ + 0(o). (31) 



By integration of this expression over the volume, the pressure mass is obtained 

3c 

When the particle moves with constant velocity the fields are 

El = qi^r (32) 

and 

B L = -vx E L , (33) 

c 

where r is the position from the centre of the particle in the laboratory frame and y 
has the same meaning in the instantaneous rest frame, 

y = r + (7 — l)(r • v)v. (34) 

The local field energy and momentum are calculated integrating, over the volume, 
the energy density = (El + BD/i^Sir) and the momentum density c~ 2 Sl, = 
cE L x Bl/(47t). The results are 

[/ e7 (l + i(-) 2 ) + 0(a) (35) 



and 



Ui . 

3 c 

Pl = \c~ 2 U cl v + 0(a). (36) 

If we add these results to the energy (mo + mp(l — 7 _2 ))c 2 7 and momentum 
(m + mp))!) of the particle we obtain the expected results for the dressed particle, 
(m c 2 + [4)7 and (m + U c /c 2 )^v. 

Why does the momentum of the field have the famous 4/3 factor? Exactly for the 
same reason why the momentum of the particle is not the bare one: because the field 
is not free and therefore it has stress. In the particle case, the stress tensor P^ v is well 
separated from the energy-momentum tensor '. The volume integrals of the time 
components of the energy-momentum tensor transform as a four-vector. The volume 
integrals of the time components of the stress tensor, that is the energy and momentum 
contributions of stress, do not transform as a four-vector. In the field case it is not 
possible in general to split the tensor T^ u into an energy-momentum tensor plus a stress 
tensor. So the volume integrals of do not transform as a four-vector. That is, the 
total energy and momentum of the field do not form a four-vector. There is nothing 
wrong with this. Relativity only requires that the total energy and momentum of a 
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confined system should form a four-vector. Indeed the effect of the stress of the field is 
cancelled by the effect of the stress of the matter. In the particular case of a single rigid 
charged body, moving with constant velocity, it is possible to obtain separated stress 
and energy-momentum tensors of the field. In the rest frame the T 00 is the energy- 
momentum tensor while the stress is T l K All the above discussion is consistent with the 
results of Singal [TTj . 

The stress tensor P^ u given by (J7|) and ()30|) is identical to the term that Schwinger 
[TU] found that had to be added to T^ u in order to obtain a four-momentum with 
electromagnetic mass U e /c 2 . We see that no ad hoc assumptions are needed, the term 
is just the pressure inside the particle. 

In order to completely solve the 4/3 problem we have to prove that (|T9*j) is verified. 
The electrostatic force density is 

/« = q 2 i 2 g(y)^fr. (37) 

It is obvious that / es does not contribute to F em because / d 3 x / es = 0, but it is the 
leading contribution to the power unbalance. In fact we are assuming that / ex does not 
stress the particle, so / d 3 x / ex • 5v = 0. On the other hand the electric component of 
the field reaction force density ff r is proportional to the acceleration, therefore it gives 
a higher order contribution to the power integral. Therefore 

Jd 3 xf-5v = J d 3 x / cs • 8v + 0(a 2 ). (38) 

To calculate this integral we note that 
dr 

Sv = ^- (39) 
dt v ; 

and that r 2 = y 2 + (^y~ 2 — l)(y ■ v) 2 . 

d i xf „.6 v= ifjd> xg (y)9M^ (40) 

= q ~Y J dS vsM^r^[(T 2 - 1)(» ■ *) 2 ] (41) 
= il [dy ^QM^Ll (42 ) 

6 J dy y dt 

as expected from (fT5j). 

5. Conservation laws of the fields 

In this section we will consider the conservation laws of energy and momentum of the 
fields. If T^ v is the momentum-energy-stress tensor of the electromagnetic field and fjf m 
is the force density acting on the particle then the conservation law is 

V M T^ = -r cm . (44) 
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As the electromagnetic field is the sum of local and radiation terms, the tensor, 
which is quadratic, contains cross terms that are products of local and radiation fields, 
rptxv _ j,llu _|_ rpfiu _|_ rp^u ^ j n t roc iucing energy densities, Poynting vectors and Maxwell 
tensors, flU)) is equivalent to 
d 

—{u h + u c + u R )+V- (S L + S C + S R ) = -f cm ■ (v + 5v) (45) 

and 

-J-(S L + 5 C + S R ) + V • (T L + T c + T R ) = -/ em . (46) 
c^ot 

Now we will integrate both equations over the whole space. One has to be careful 
because decays as r~ 2 . So we take a volume V, surrounding the particle, so large 
that T£ v « and Tq & outside V. In these conditions the radiated power and the 
rate of radiated momentum defined as follows are independent of V. 

Prad = ^ / d 3 x« R + / dA-S R , (47) 
at Jv Jdv 



We get 



Grad = d 6 xS R + dA- T R . (48) 

C Clt JV J dV 

^(U L + U c ) + P rad = ~F cm ■ v - (49) 
at 6 dt 

^-(P L + P c ) + G rad = -P em . (50) 
dt 

Here we have used (I4UJ1 and the total cross energy and momentum defined as 

[/ c = / d 3 xu c (51) 



and 



and 

Pc = - 



l -jd*xS c . (52) 

The dressed particle includes, in addition to the matter momentum and the pressure 
momentum, the momentum of the local fields. The reaction of these is Pl = — Pl- This 
reaction should be subtracted from the self-force P em in order to obtain the force that 
the dressed particle feels. We will call this last force radiation reaction, 

dPL 

Pr = Pern — Pl = Pern + 

Actually, F R contains, not only the reaction of the radiated fields but also the 
reaction of the cross terms, 

Pr = -Grad - ^jjp (54) 

The power of Pr is obtained from (f4T?|) using (jHSJ), and the identity (fTHJ) . It also 
contains a contribution from the cross terms, 

F R -v = -P iad -^. (55) 
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Finally, the motion equation and the power equation of the dressed particle are 
readily obtained from (JTljl and (jl7j) . with dressed mass m = m + U e /c 2 

mc 2 ^ = (F cx + F R )-v. (57) 
6. Radiation reaction 

The calculation of the self-force is much simpler if it is done in the reference frame 
of instantaneous rest. To begin with, there are no magnetic field contributions so we 
have only to calculate the electric field. It is very convenient in this case to split 
the electric field into electrostatic and induced components, E = E e + E h with the 
following properties: V • E e — 47ip, V x E e = and V • E\ = 0. Note that E e is 
the longitudinal and E\ the transverse field. With these definitions, in the rest frame, 
/ es = pE c and /f r = pE L . As we have previously explained in the integrated self-force 
-F cm the electrostatic contribution cancels out. The easiest way to determine these fields 
is to use the Coulomb gauge, V • A = 0, because in this case the potential <fi yields E c 
and A yields Ei, 

E e = -V0 (58) 



and 



E, = - l - S 4. (59) 



c dt 

The electrostatic field is 



p{y,t) 



E c (x, t) = I d 3 y P^ (x - y) (60) 
\x — y\ 



and A can be obtained from the equation 

1 d 2 A 2 A 4tt . , 1 dE c 



_V 2 A = — j + -^=r. (61) 



c 2 dt 2 c c dt 

The solution is 

Aix , t) = if A s y fMl + i / d 3 i W 

cJ \x — y\ Aire J \x — y\ ot 
where t' is the retarded time, t' = t — \x — y\/c. The induced field is therefore 

*<-*>- -hl^ 9 -^ - ^I^ 9 ^ 

For the quasi-rigid motion the charge density is p(x,t) = q r yg{y) and the current 
density is j = p(v + 5v), where y = r + (7 — l)(r • v)v, y = \y\ and r = x — Xo(t). 
With these definitions we have in the rest frame 

d 2 p(x,t) dq A . . ^ . 9n . aS 

= ~ q I/ ' ait) + ° (a } (64) 
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i cp. f \ 

= qg{r)a(t) + 0{a 2 ). (65) 



dt 

If one now replaces the expression for p into the electrostatic field equation one gets 

d 2 E c {x,t) _ d_ f d z dg V ■ a(t) ^ 
dt 2 dx J dy \x — y\ 

The integral in (jHEJ) can be evaluated, 



* (67) 
ay | x — y\ x z 



here Q(x) is defined in (I20j) and we have used the fact that 
1 fjn y min(x,y) . 

— / dSL = —x. (68) 

An J \x — y\ 3max(x, y) z 

The second derivative of the electrostatic field is obtained by taking the gradient 
of (|EI|), then 



d EJxA) r , * * * Q(x) , . , 

V ' - - 4Tig(x)xx + -^(I-3xx) ■ a(t). (69) 

x s J 



-9 



F m (*) = ~4 f[d s xd> V &^a{t-\x-v\/c). (72) 
6c z JJ \x — y\ 



dt 2 

Here I is the identity tensor. This expression should be replaced into (|53*|) and then 
integrated again in order to obtain the self-force 

F cm (t) = qJd 3 xg(x)E i (x,t). (70) 

In this process the induced field is spherically averaged. Because of (J65|) and as 

— / dQ x xx = -I, (71) 

one finds that the contribution of the second term of (J63|) is exactly —1/3 of that of the 
first term. Finally we obtain the self-force formula in the rest frame 

2( f rr J 3- J 3.. g( a Qg(y) 

\x - y\ 

As was expected it is an average of the retarded accelerations, therefore no violation of 
causality is due to this formula. By using the Coulomb gauge we have obtained F em 
in closed form instead of the troublesome power expansion one gets when the Lorenz 
gauge is used. The sum of that power expansion appears in the textbook of Jackson 
jHj, and indeed our equation (J72J) can be deduced from equations (16.28-16.30) of that 
book. 

When the acceleration is constant the self-force reduces to the reaction of the local 
fields, 

= (73) 
One has to subtract this contribution in order to obtain the radiation reaction force 

= 1-2 II ^ d "y 9 -^l("(t) - a(t -\x- y\/c)). (74) 
Sc z JJ \x — y\ 
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If one expands the retarded acceleration in terms of \x — y\ the LAD result is 
recovered as the first term in the expansion, 

F R =l^d(t) + 0(R). (75) 

If one uses (172|) to calculate the self-force of a charged spherical shell of radius R 
the Sommerfeld-Page-Caldirola-Yaghjian result is obtained, 

F cm (t) = ~-^(v(t) - v(t - 2R/c)). (76) 

7. Integral equation and point particle equation 

In the rest frame the equation of motion is 

ma = F cx + F R . (77) 

This is, of course, the non-relativistic equation. The condition a <C c 2 /(2R) implicit in 
the quasi-rigid approximation implies that « < c for times of the order of the transit 
time (i. e. the time that the light takes to traverse the particle) and therefore the non- 
relativistic equation can be used. Our purpose is to solve (|77j) in order to express the 
radiation reaction as a functional of the external force, instead of a functional of the 
acceleration. Equation (J73J) can be written as 

/•oo 

F R (t)= dt'k{t'){a{t)-a{t-t')) (78) 

JO 

where the kernel k is given by 

k(t) = ^- ffd 3 xd 3 y 9 ^ )9 ^ 5(t-\x-y\/c). (79) 



3c 2 J J \x — y 

The solution is obtained by means of a Fourier transform 



K(u) = Tdt k(t)e^ = % ffd 3 x d 3 y ^iv) ^-v\/^ (80) 
o 3c 2 J J \x — y\ 



The mass of the local fields equals K(0), 



2 1 2 /Y j3 i3„ g(x)g(v) 4 



while the derivative of K(oS) at uj = gives the constant of the LAD equation, 
.dK f°° 2r/ 2 /•/• ,,, , ;! 2g 2 



du 



o 



J o °°dttk(t) = -L J J d 3 xd 3 yg(x)g(y) = (82) 



From ()81j) and (|H2*|) one obtains an expression for the average delay time r 



1 



g(x)g(y) \ 1 



m j a ^^{cJJ^W^) ■ (83) 

The acceleration, obtained from (|77|). is then 

/oo 
dt'G(t')F cx (t-t'), (84) 
-oo 
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where 

1 roo p—'iut 

G(t) = — duo — — -. (85) 

w 2tt J-oo m- K{0) + K{uj) v ; 

Causality requires that G(t) = for t < 0. That implies that (m — K(0) + K{uj))~ 1 
must be regular for lm(u) > 0. K(u) is regular for lm(u) > 0, so the only possible 
singularities are poles when m — ^(0) + K(oS) = 0. For rj > K{xrj) is real and 
lim^oo K(irj) = 0. When g(r) > 0, K(ii]) decreases monotonically from K(0) to 
zero. Even when there are regions in which g(r) is negative K(ir]) > 0. Therefore 
(m — K(0) + K^uj))" 1 is regular in the upper half of the complex plane if, and only if, 
m > K(0). That is, it must be m > ^U c /c 2 , or equivalently m > ^U e /c 2 . Note 
that this is the same condition discussed in section 3. If the condition is fulfilled 
G(t) ~ 6{t) exp(— t/ra) and tq ~ t. If the condition is not satisfied, causality is 
violated and run-away solutions appear. Let us look in more detail how this happens. 
Consider a case in which the external force was constant for t < 0, and then suddenly 
vanishes at t — 0. The self-force F em given by (|72|) is opposite to the acceleration and 
remains the same at the transition, F em (0 + ) = F em (0~) = —K(0)a(0~). The motion is 
determined by (JHJ). For t > the acceleration equals F em divided by the effective mass 
of matter m + mp = m — K(0), in particular «(0 + ) = — (m — K(0))^ 1 K(0)a(0^). If 
itlq + mp > the new acceleration is opposite to that of negative times and, as time 
elapses, the accelerations for t > subtract from the previous ones in ([72)1: therefore 
the self-force is reduced and the acceleration vanishes exponentially. On the contrary if 
< m < K(0), the new acceleration not only has the same sense of the previous one but 
it is actually larger than it. As time elapses the self- force increases and the acceleration 
diverges exponentially. 

The mass condition was previously demonstrated for the spherical charged shell by 
Moniz and Sharp Although their quantum-mechanical results have been disputed 
[2*U] . their classical calculation of the run-away solutions is correct, apart from the fact 
that they do not include the pressure contribution to the mechanical mass. Here we have 
shown that the mass condition is valid for any charge distribution and have clarified its 
relation with the inertia of stress. 

The equation ()84jl can be re-written as 



a{t) 



roo roo 

/ dt'G(t') F m (t)+ / dt'G(t')(F cx (t-t')-F cx (t)). (86) 
Jo J Jo 



The pre-factor of the force in the first term of ()86)1 happens to be m 1 



oc 



dtGCt) = I du -4v \ = — . (87) 

The radiation reaction is then 

/■CO 

F R = m dt' G(t')(F cx (t - t') - F cx (t)). (88) 
Jo 

Finally, if the external force varies slowly enough, more precisely if 

^ « r- 1 (89) 
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an approximated instantaneous radiation reaction is obtained by taking the first term 
of the expansion in t' of the retarded force in 

„ 2 q 2 ■ 



3 c A m 

To obtain this we have used the fact that 

1 d8(u) 



(90) 



roc poo 

/ dttG(t) = i / du 

JO J-oo 



m-K(0)+K(uj) duo 



i dK 



m 2 duj 



2 q 



2 



3 c 3 m 2 



(91) 



The equation (|90j) was proposed by Rohrlich J7| as the correct one for point 
particles. We see that as long as the mass condition m > ^U e /c 2 is fulfilled the radiation 



reaction can be expressed by (j88|) which conforms to causality and does not produce 
run-away solutions. The approximated expression of (J90|) does produce violations of 
causality when Fr ~ F ex . When this happens the condition of (fHH|) is not satisfied, so 
the approximation is not valid. 

A particle can be considered a point particle if its radius is negligible when compared 
with the wave-length of the radiated field. This condition is equivalent to ([89)1 . So (j9*Uj) 
can indeed be considered the correct one for point particles. It is worth noting that 
applied forces with discontinuities, that violate (|89p. are very common. However, usually 
such discontinuities are not real, but the result of some approximation. In most cases 
the formula of (j97Ijl behaves properly if, instead of discontinuities, realistic smoothed 
steps are used. 



8. Radiation reaction four-force and radiated power 

The radiation reaction of the physical point charge with v = is given by equation (j90j) . 
In this section we find this force for a particle moving with arbitrary speed. 

A four-force that conserves mass is related to the force by K % = 'jF 1 and to 
the power by K° = • v/c. This implies that it is orthogonal to the four- velocity 
K^u^ = 0. In the rest frame K l = F % and K° = 0. 

The motion equations of the dressed mass, and (jKTjl . show that the radiation 
reaction is a mass-conserving force. Then, the four-force of the radiation reaction of the 
physical point charge is in the rest frame 

K*=- i ex (92) 

R 3c 3 m dt V ; 

and 

K° R = 0. (93) 

Therefore, if we call the external four-force and r the proper time, the radiation 
reaction four-force is in an arbitrary reference frame given by 
2 q 2 



... - j 2 vdK^ 1 AK V 

Kr 



3 (?m 



-z — : — U, 



dr c 2 dr 



(94) 
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The four-scalar in parenthesis is 



a. 



From (|94|) the force and power expressions are obtained for any speed. 



R 



2_q 2 r d 7 F ( 
3? 



m 



r 



dt 



a v 



and 



jp 2 q 2 r d 7 F { 



3 c 3 m 



df 



7^ 



(95) 



(96) 



(97) 



From these equations, and equations (|54j) and (|55|) one obtains the momentum and 
energy of the cross terms, the radiated power, and rate of radiated momentum, 

2 q 2 



Pc = 

U c = 
p 

1 rad 



3 c 3 m 

1 A 2 

3 c 3 m 

2 q 2 



7^ex, 
7^ex • v 



P c -v, 



YF C , 



3 c 3 m 
2 g 2 7 3 



a, 



0)« = ^rad^- 



(98) 
(99) 
(100) 
(101) 



3 c 3 m c 2 

Equation p00|) reduces to the Larmor result ^2] m the non-relativistic limit only 
if i^x • a = 0. The Larmor result is valid for a mathematical point charge (R = 0), not 
for an extended particle. 



9. Conclusion 

We have completely solved the self-force of the extended quasi-rigid particle. We 
have shown that, provided m > |f/ e /c 2 , it is a perfectly consistent classical system 
conforming to causality and conservation of energy and momentum. In order to obtain 
the correct equations it is essential to use the relativistic conservation equations of energy 
and momentum for the continuous medium, even if v c; only so the inertia of pressure 
is properly taken into account. The mathematical point charge (R = 0) is inconsistent 
with classical electrodynamics. For a given dressed mass m and charge q, the minimum 
radius that any particle can have is (2g 2 )/(3mc 2 ). Nevertheless if the radius of a particle 
is much smaller than any other distance in the problem, the radius becomes irrelevant 
and the particle can be properly considered as a point charge from the physical point 
of view. 

A comment about QED. A correct quantum theory of a system, should give in the 
classical limit a proper classical system. As the mathematical point charge is inconsistent 
in classical electrodynamics, it follows that the proper quantum theory of a structureless 
particle (such as the electron) should have as classical limit an extended classical particle. 
At first sight one might think that this were not impossible as the single particle picture 
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breaks down at the Compton length scale, which is a -1 (~ 137) times larger than 
the classical radius of the particle. The relativistic quantum non-localities produce an 
effective particle extension. Nevertheless the QED calculation for spin | of Low [21] 
shows that these effects are not enough to eliminate the run-away solutions, but that 
the size bound is reduced by a factor of exp(— a -1 ). 
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